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Abstract 



A recently proposed parametrization for the deep inelastic diffractive cross section is 
used to describe the HI 94 data. We find two possible solutions, and we discuss them is 
some detail. 



1. Diffractive events at HERA fl], § are generally attributed to the exchange of the 
Pomeron between the virtual photon and the proton. After the first observation of these events, 
attempts of a theoretical description || where based upon the idea of Ingelman and Schlein 
0, and they viewed diffraction in deep inelastic scattering as "deep inelastic scattering on a 
Pomeron". As a striking result along these lines, HI |J reported on a "hard gluon component 
of the Pomeron" (where "hard" refers to the behavior of the structure function near (3=1). 
However, with increasing statistics both ZEUS and HI [jl], [!| found evidence that in DIS 
diffraction the Pomeron intercept tends to be higher than the soft hadronic Pomeron. This 
observation suggested that the diffractive cross section at HERA contains both a soft and a 
hard component, and the latter one is stronger than initially expected. 

A theoretical description of the inclusive diffractive cross section, therefore, should contain 
both a nonperturbative Pomeron part and a perturbative contribution. As a possible strategy, 
one could start from that part of the final state where pQCD can be used safely, and then 
attempt to find a smooth extrapolation into the nonperturbative Regge region. As to the 
Pomeron, in the perturbative region it is most easily modelled by a two gluon-system, the gluon 
structure function of the proton. In the nonperturbative region, this model of a "hard Pomeron" 
should smoothly connect with the soft hadronic Pomeron. Recently || a parametrization of 
the diffractive cross section has been proposed which follows this line of arguments, and it has 
been applied to both ZEUS and HI data. For the latter the kinematic region was restricted to 
those values of (3, Q 2 , and xjp for which also ZEUS data points exist. As a striking result, it 
was found that the ZEUS data allow for a unique and rather succesful description whereas for 
the HI data two different fits of this model were found. 

In the present paper we apply the parametrization of || to the complete HI 94 data of 
the diffractive cross section. We confirm the existence of two different solutions and we study 
their properties in some detail. 

Let us first briefly recapitulate the parametrization. It consists of 3 terms each of which 
stands for a particular diffractive final state: 

F^^A^Y'pil-P) (1) 

^ (3 ,,„ = B gy as ^ ( | + 1) j (1 _ /J)7 (2) 
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where 

n 2 , 4 = + n\j,n 

The first term describes the diffractive production of a qq pair from a transversely polarized 
photon, the second one the production of a diffractive qqg system, and the third one the 
production of a qq component from a longitudinally polarized photon. The ansatz for the 
/3-dependence is motivated by rather general features of QCD-parton model calculations: At 
small (3 (large diffractive masses Mx) the spin 1/2 (quark) exchange in the qq production (1) 
leads to a behavior ~ /3, whereas the spin 1 (gluon) exchange in (2) corresponds to j3°. For 
large (3 (small diffractive masses) in (1) and (3) calculations show || |7j that pertubative QCD 
becomes reliable, and they lead to 1 — /3 and (1 — (3)° in (1) and (3), resp. For the qqg term 
in (2) the situation is slightly more complicated: a QCD calculation [§] leads to a (1 — /3) 3 
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behavior. But in view of the previous HI result which suggests a hard gluon distribution inside 
the Pomeron and, therefore, also a large gluon contribution in F 2 near (3 = 1, we will leave 
the exponent 7 as a free parameter. Comparing the behavior in (3 near 1 of the three terms we 
notice that the longitudinal term (3) dominates. 

Turning to the Q 2 dependence, the first two terms are leading twist whereas the last one 
belongs to twist 4. Combining this with the (3 dependence discussed before we see that the 
longitudinal term is important only near (3 = 1. The lnQ 2 -terms follow from QCD-calcutions, 
and they indicate the beginning of the QCD Q 2 -evolution. In the transverse case ((1) and 
(2)) the evolution starts only after the radiation of at least one gluon, whereas for the higher 
twist longitudinal case the first logarithm appears already in the quark loop. The longitudinal 
qqg production contributes to leading twist, but it does not have the logQ 2 enhancement and 
therefore is small compared to the transverse term (2). It should be clear that the parametriza- 
tion (1) - (3) does not yet aim at describing the Q 2 evolution over a large Q 2 region. At this 
early stage it is designed only for the small-x HERA region where Q 2 stays below 100 GeV 2 . 
The use in the large-Q 2 region will definitely require some refinement of the parametrization. 
Finally, the dependence on xjp cannot be obtained from perturbative QCD and therefore is left 
free. For the two transverse terms one expects approximately the same xjp-behavior, not too 
far away from the soft Pomeron. For the longitudinal term (3), on the other hand, it has been 
shown [[?] that the xjp dependence should be approximately the same as for the square of the 
gluon structure function. In the fit both 112 and are allowed to have have a linear rise in 
lnlnQ 2 . 

In contrast to the ZEUS diffraction data HI data extend into the region of not so large 
rapidity gaps where, in the Regge language, also the exchange of secondary reggeons has to be 
taken into account. Within the spirit of the parametrization proposed in || such a term is 
very natural: within pQCD it corresponds to the exchange of a qq system. Such a term has 
recently been studied in 0. In our present fit we follow the HI procedure of [JTJ] and use the 
following ansatz: 

F 2 D ^ IV ( XlP ,(3,Q 2 ) = Nf R/p (x M )F 2 M ((3,Q 2 ) (5) 

where the reggeon flux is taken to follow a Regge behaviour with a linear trajectory a R {t) = 
a K(0) + ct R t such that: 

f-tmin e B R l 

fjR/pM = 2a M (t)-l dt ( 6 ) 

"'tout Xjp 

where \t m i n \ is the minimum kinematically allowed value of \t\ and t cut = — 1 GeV 2 is the limit 
of the measurement and the values of Bm and a K are fixed with hadron-hadron data jl]. F 2 R 
is assumed to be proportional to the pion structure function |10[ []. No interference is assumed 
between this reggeon exchange and the Pomeron exchange of (1) - (3). The presence of (5) is 
a new element in our fit. 

In our fit the main parameters are the normalization constants A, B, C, N, and the 
exponent 7 of the /3-dependence near (3 = 1. Moreover, we allow the xjp exponents to vary, i.e. 
we have the five parameters n 2 , n\, n®, n\, and the intercept an(0) of the secondary. Finally, 
we also leave free parameter. 

2. Let us now turn to the results of our fit (table 1). Most important, the data allow two 

*It was checked that assuming another pion structure function JTT| for the reggeon does not modify the 
results of the fits. 
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Table 1: Parameters obtained for the "perturbative gluon" (on the left) and "hard gluon" (on 
the right) solutions of the fit. The fits are performed with statistical and systematical errors 
added in quadrature. The two last fits are preformed imposing the qqg (B=0) or the longitudinal 
qq components to be zero (C=0). If the fits are performed with statistical errors only, the x 2 
values for the perturbative gluon and for the hard gluon solutions are 236.0 (1.10/dof) and 
269.1 (1.25/dof), resp. 

different possible fits of the 1994 HI data[] The parameter which distinguishes the two solutions 
is 7: in the first case we find 7 = 8.24, in the second case 7 = 0.27. Both solutions have an 
acceptable x 2 value, with a slight preference for the first one (high gamma). If we perform the 
fit with statistical errors only, the x 2 values per degree of freedom of the two fits are 1.10/df 
and 1.25/df, resp. An overall picture of the two solutions is given in Figures 1 and 2. 

The difference between the two solutions is most clearly illustrated in Figs. 3 and 4 where 
we plot the /^-dependence of the different pieces (1) - (4). Data points lie mostly in the central 
and the rhs columns (xjp = 0.001 and xjp = 0.01), but for illustration we also extrapolate 
down to X]p = 0.0001. The main difference lies in the behavior of the qqg component: for 
the first solution (7 = 8.24) it contributes only at small (3 {(3 < 0.3), whereas for the second 
one (7 = 0.27) it extends up to (3 near one. Qualitatively this second solution is close to the 
previous HI solution with the hard gluon: the fact that the gluon distribution obtained from 
the HI QCD analysis peaks near 1 leads to a large gluon contribution to F 2 D near (3 = 1. In 
the following this solution to our fit will therefore be referred to as the "hard gluon" solution. 
The first solution, on the other hand, one is closer to the perturbative prediction which leads 
to 7 = 3. This solution will be called "perturbative gluon" solution in the following. For 
both solutions one sees that the transverse qq term has its maximum at medium (3, and the 
longitudinal one contributes only near (3 = 1. For the perturbative gluon solution the transverse 
qq is somewhat bigger than in the hard gluon case. 

An important feature of both solutions is that none of the four components (1) - (4) comes 
out to be particularly small. To illustrate the importance of (2) and (3) we have repeated the fits 
(table 1) with putting either B or C equal to zero: in both cases the x 2 increases significantly. 
Let us notice that only one solution remains if B = is required (both solutions in Table 
1 are identical). On the other hand, imposing C = in the hard gluon solution implies a 
very low value of 7 (7 = 0.001), and so an important qqg term is present even at high (3 and 
compensates the smallness of the qq term in this region. The situation changes if we impose 

^To be more precise we mention that there are two more solutions: they differ from those in table 1 by their 
value of xq. For the perturbative gluon solution we find xq = 1.9 10~ 4 , for the hard gluon solution xo = 4.9 10~ 5 . 
Their x 2 values are substantially higher: \ 2 = 242.2 and % 2 = 223.8 for the perturbative and for the hard gluon 
solution, resp. Therefore we will not discuss them in our paper. 



HI 1994 Data 




Figure 1: Result of the "perturbative gluon" fit. The triangles are the HI 1994 data. In full 
line is displayed the full result of the fit, and in dashed line, the reggeon component only which 
is important at low (3 and high xjp. 



HI 1994 Data 




Figure 2: Result of the "hard gluon" fit. The triangles are the HI 1994 data. In full line 
is displayed the full result of the fit, and in dashed line, the reggeon component only which 
is important at low (3 and high xp. This hard gluon fit is hardly distinguishable for the 
perturbative gluon one. 
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Table 2: Parameters obtained for the perturbative gluon solution of the fit after different cuts 
on (3 in the data. 

the same kinematic cuts as used in the QCD fit of [0]. Namely, by excluding the data points at 
high y (y >0.45) and at low Mx (Mx <2 GeV) this fit was restricted to the kinematic region 
where higher twist corrections are expected to be small. Also, only (3 and Q 2 bins with more 
than three measured points were considered. Imposing the same cuts we have fitted our model 
to the data: we still find the two solutions. However, as expected, the longitudinal qq is now 
poorly constrained, and forcing the C parameter to be in the fit (excluding the higher twist 
contributions) does not change the quality of the fit very much (for the perturbative gluon 
solution the x 2 increases from 130 to 133, whereas for the hard gluon it remains unchanged 
(x 2 =144)). The hard gluon solution now is even closer to the previous HI solution in which 
no qq term was present. 

The importance of the secondary reggeon term can be seen from Figs.l and 2: in both 
cases its contribution becomes substantial only at low (3 and large x$>. 

Scaling violations are illustrated in Figs. 5 and 6. The Q 2 dependence of the longitudinal 
qq term (3) near (3 = 1 is easily understood: it decreases with increasing Q 2 , but the decrease 
is much slower than 1/Q 2 : this is due to the (logQ 2 ) 2 enhancement in (3). As to the other 
two terms (1) and (2), there is an additional Q 2 -dependence coming from the exponent n^- in 
both fits n\ is not small, and in combination with the rather large values of the ratio Xq/xjp it 
leads to a substantial increase with Q 2 . In Figs. 5 and 6 this explains why also the transverse 
qq component rises with Q 2 . For the qqg component, which according to (2) has a logarithmic 
growth in Q 2 , we see a difference between the two solutions: in the perturbative gluon solution 
it contributes only at small (3. Consequently, the rise in Q 2 of is stronger at small (3 than at 
larger (3 (e.g. (3 = 0.65). This effect is absent in the hard gluon solution, since here the gluon 
is present at all (3. 

In order to gain futher insight into the two solutions, we have also performed fits in a 
limited range in (3 (or Mx)- The parameters obtained for the perturbative gluon fit with cuts 
on (3 in the data are given in Table 2. We note that the parameters are quite stable except 
for 7 which has the tendency to increase. However, the error on 7 is large, because we cut 
explicitely the region at low (3 (or large Mx) where the qqg is significant. The hard gluon fit, 
on the other hand, is much less stable: already with the lowest cut at /3 > 0.1 this solution 
disappears, and only the perturbative gluon solution remains. Cuts on Mx in the data lead to 
the same results: excluding the high Mx region, the perturbative gluon solution turns out to 
be stable, whereas a cut Mx < 10 is already enough to eliminate the hard gluon solution. 

3. Next we use our ansatz (1) - (4) for the ZEUS 94 data @. In agreement with [|| 
we find only one solution, and the 7-value turns out to be identical: 7 ~ 4.3. Comparing 
this solution with the perturbative gluon solution of our HI fit, we find that they are very 
similar, both qualitatively and quantitavely. For instance, if in the HI fit we fix 7 at the ZEUS 




Figure 3: Display of the different components for the "perturbative gluon" fit. The transverse 
qq term belongs to the dashed line. The qqg term (dashed-dotted line) is important at low f3, 
and the longitudinal qq term (dotted line) dominates at high (3. 




Figure 4: Display of the different components for the "hard gluon" fit. The transverse qq term 
belongs to the dashed line. The qqg term (dashed-dotted line) is not much /9-dependent, and 
the longitudinal qq term (dotted line) dominates at high j3. 
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Figure 5: Scaling violations in (3 bins for a fixed value of xp = 0.003 for the perturbative 
gluon (on the left) and hard gluon (on the right) solutions of the fit. The transverse qq, qqg, 
and longitudinal qq components are respectively in dashed, dashed-dotted and dotted lines, 
respectively. The full line is the sum of all components, including the secondary reggeon. 



value 7 = 4.3, we get a x 2 °f 197.8 which is not much bigger than the \ 2 °f the fit with all 
parameters left free (x 2 = 186.7). As an important result of our ZEUS fit, we note that the 
fit is not improved by adding the secondary reggeon: the reggeon normalisation is found to be 
compatible with zero. This shows that ZEUS data do not require any secondary contribution. 

To obtain a better understanding of the differences between the fits to the HI data and 
to the ZEUS data , we have tried to compare the result of the HI fit (the perturbative gluon 
solution) with the ZEUS data. The result is shown in Fig.[7]. The triangles represent the recently 
published ZEUS data, and the squares the HI data. To be able to compare directly HI and 
ZEUS data, we had to shift the HI data to the Q 2 and (3 ZEUS bins, using the perturbative 
gluon solution of the fit. Let us notice that this comparison does not depend much on the 
way we perform the extrapolation, as the Q 2 and (3 bins from HI and ZEUS are rather close. 
Further more, an extrapolation using a completely different fit based on the dipole model give 
the same interpolated values fl2fl . The full curve is then the perturbative gluon fit result with 
HI data. The hard gluon solution restricted to the ZEUS bins would be indistinguishable on 
this plot. The dashed line denotes the result of our fit to the ZEUS data. We notice that the 
fits differ precisely in that region where one sees differences in the data, more precisely in the 
bins (3 =0.2, Q 2 =8, and, j3 =0.7, Q 2 =60 [12|. The differences in the fit parameters are clearly 



due to differences in the data points. More precise data in these regions will be thus very useful. 
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Figure 6: /3-dependence of F 2 D in bins of Q 2 at a fixed value of xjp = 0.003 for the perturbative 
gluon (on the left) and hard gluon (on the right) solutions of the fit. The dashed, dotted and 
dashed-dotted lines represent the transverse qq, the longitudinal qq (important only at high 0), 
and the qqg terms. The full line is the sum of all components, including the secondary reggeon. 



4. So far we have not discussed our results for the exponents n<i and n±. Since the errors 
for these parameters are not small, one may feel that at this stage we should not attribute 
too much weight to them. Nevertheless, in the following we describe a few fits in which we 
impose constraints on n<i and/or n 4 . First we try to identify the xjp dependence in (1) and (2) 
with the soft Pomeron, leaving as a free parameter. We put 712 = 1.12 (this differs from 
the soft Pomeron value 2ajp(0) — 1 = 1.16 by a small amount which is due to the integral 
over the momentum transfer t) and n\ = 0, i.e. we exclude any Q 2 variation. The results 
for both solutions are shown in table 3. The obtained x 2 are quite bad: 254.2 (291.2) to be 
compared with 206.8 (184.6) in table 1 for the perturbative (hard) gluon solutions. This makes 
it unlikely that the diffractive structure function data could be described by the soft pomeron. 
In a variant of this test we have also treated n° as a free parameter (still keeping n\ — 0). For 
the perturbative gluon solution we find a reasonable fit (x 2 = 221.9) for n\ = 1.38, and for the 
hard gluon solution the fit gives = 1.23 with x 2 = 239.6. 

Next we test the possibility of identifying the xjp dependence of the longitudinal term 
with the gluon density in the proton. We take the gluon parametrisation from [ll which is 
known to give a good description of the proton structure function data measured at HERA and 
fit it to the form x~ a+1 . We find 



xG ~ x 



-a+l 



-0.18-0.15 ln(ln(Q 2 )+l) 



(7) 



i n 




Figure 7: Comparison between HI (squares) and ZEUS data (triangles). The full line represents 
the HI perturbative gluon fit, and the dashed line our ZEUS fit. We notice that the differences 
in the fits are precisely in the region where the data differ. 
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where the scale of Q 2 is 1 GeV 2 . This suggests to parametrize the exponent in the following 

way: 



n 4 = 2a - 1 = 1.45 + 0.171n ( hi(^) + 1 



q 2 



The result of both the perturbative gluon solution and the hard gluon solution are given in 
Table 3 (middle columns). We note that the quality of the fits is about the same as those 
obtained in Table 1. The transverse terms prefer a soft Pomeron with a rather strong Q 2 
dependence of the exponent n 2 . We also observe that, for the hard gluon solution, the value of 
the parameter xq is quite different from the fit in table 1 (cf. the footnote on p. 2). We conclude 
that both solutions are consistent with the theoretical expectation that the longitudinal part of 
our parametrization is hard, i.e. its energy dependence can be described by the gluon structure 
function. 

We have also tested the possibility of identifying the xjp dependence of the longitudinal 
term with the structure function F 2 of the proton. It has been shown |13| that the F 2 slope in 
x can be parametrised as follows: 

where the scale of Q 2 is 1 GeV 2 . The exponent can then be expected to vary as 



Q 2 
Ql 



n 4 = 2a- 1 = 1.31 + 0.18 log 10 ^ ( 10 ) 



In our fit we put the scale Q\ = lGeV 2 . The results are shown in the righmost column of Table 
3. Again, the transverse terms in our parametrization want to have a soft Pomeron with a 
rather large Q 2 dependence. The quality of the fits is about the same as that of the first fits in 
Table 1. Changing the scale to Q\ = 4GeV 2 or treating it as a free parameter does not change 
the result very much. Alltogether, the parameters of these fits are very similar to the previous 
one with the gluon structure function. 

5. We conclude with a few general comments. Within the parametrization suggested in 
II have tried to clarify the differences between the two possible solutions for the HI 94 data, 
and their relation to the previously reported HI solution with the hard gluon. One of the two 
solutions (named 'perturbative gluon') is closer to what one expects from a straightforward 
extrapolation of perturbative QCD calculations, whereas the other one (called 'hard gluon') 
is closer to the previous HI result. The main difference lies in the f3 dependence of the qqg 
contribution. Apart from a small difference in y 2 (in favor of the perturbative gluon solution), 
we also found some evidence that, when we impose kinematic cuts on the data the hard gluon 
solution is slightly less stable than the perturbative gluon solution. 

As we have discussed in the beginning of this letter the parametrization for F 2 D is based 
upon an attempt to interpolate between the QCD parton model (hard region) and Regge physics 
(soft Pomeron). In this first step of testing this approach we have given much freedom to our 
parametrization. In a future step it will be interesting to see how specific final states fit into 
this parametrization and maybe used to constrain the parameters of our fit. Two obvious can- 
didates are the production of longitudinal vector mesons which should be contained in the third 
term of our parametrization, and the production of qq jets (contained in the first term) and qqg 
jets (inside the second term) from the transverse photon. The measured cross sections of these 
final states and a careful analysis of their topological features could be used to estimate the 
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0.30 ± 0.19 






XQ 


0.19 ± 0.03 


0.10 ± 0.04 


0.12 ± 0.02 


x A 


254.2 


190.9 


186.7 



7 


0.70 ± 0.11 


0.16 ± 0.11 


0.23 ± 0.13 


<*R 


0.40 ± 0.04 


0.75 ± 0.03 


0.76 ± 0.03 


N 


38.6 ± 4.8 


7.9 ± 1.5 


7.3 ± 1.3 


A 


0.086 ± 0.001 


2.01 ± 0.51 


1.16 ± 0.29 


B 


0.045 ± 0.005 


0.75 ± 0.21 


0.37 ± 0.11 


C 


0.018 ± 0.011 


1.70 ± 0.67 


1.28 ± 0.42 


< 


1.12 (fixed) 


1.04 ± 0.09 


1.05 ± 0.06 


nl 


0. (fixed) 


0.30 ± 0.07 


0.29 ± 0.05 


J 
t 


1.44 ± 0.06 






0.09 ± 0.09 






XQ 


0.41 ± 0.04 


0.014 ± 0.002 


0.022 ± 0.003 


x 2 


291.2 


210.0 


208.3 



Table 3: Parameters obtained for the perturbative gluon (on the left) and hard gluon (on the 
right) solutions of the fit. The first fit is obtained by imposing the value of the soft pomeron 
exponent for n 2 and in the second and third fit respectively, n 4 follows the gluon slope in x and 
the F 2 slope in logl/x (cf text). 

parameters of our ansatz. Moreover, their measurement can decide which of the two solutions 
discussed in this letter is to be preferred. 

Acknowledgement: We thank J.Dainton for a careful reading of our manuscript and for 
his helpful criticism. We also thank R.Peschanski for a careful reading of the manuscript. 

Note added in proof: Another fit of our model to the combined 94 and the (preliminary) 
95 low Q 2 data of HI has been presented at the Brussels DIS98 (Talk presented by T.Nicholls) 
conference. The extrapolation into the region of smaller Q 2 values works quite reasonably. 
More recently (HI Collab., Contribution to ICHEP98, Vancouver, July 98), our fit has been 
extrapolated also into the large Q 2 region (200 GeV 2 < Q 2 < 800 GeV 2 ). We feel rather hesitant 
in using our simple parametrization in this region. First, as we have said before, extending the 
Q 2 region in such a dramatic way our parametrization needs more terms and a more accurate 
way of handling the Q 2 evolution. Secondly, in the large Q 2 region region of HERA xp is not 
really small, and the influence of the secondary exchange term (4) of our parametrization starts 
to dominate. Comparison of our fit with the large-Q 2 data, therefore, does not only test the 
Q 2 evolution of our parametrisation but also of the pion structure function used in the fit. 
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